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Abstract 

Important features of low-frequency collective vibrational excitations in neutron drip line nuclei 
are studied. We emphasize that pairing anti-halo effect in the Hartree-Fock-Bogoliubov (HFB) theory 
play crucial roles to realize collective motions in loosely bound nuclei. We study the spatial properties 
of one particle - one hole (lp-lh) states with/without selfconsistent pairing correlations by solving 
simplified HF(B) equations in coordinate space. Next, by performing Skyrme-HFB plus selfconsistent 
quasiparticle random phase approximation (QRPA) we investigate the first 2 + states in neutron rich 
Ni isotopes. Three types of calculations, HFB plus QRPA, resonant BCS plus QRPA, and RPA are 
compared. 

1 Introduction 

The study of the properties of nuclei far from the /?-st ability line is one of the most active and challenging 
issues in nuclear structure physics. Especially interesting phenomena concern the loosely bound neutron 
systems close to the neutron drip line. Low-frequency collective vibrational excitations in neutron drip 
line nuclei are one of the most interesting subjects. Naively we expect low-frequency vibrational modes 
associated with low-density neutron matter (neutron skin and halo). A famous one is soft dipole mode 
012). However it is a fundamental question whether such collective modes can really appear or not. 

Vibrational excitations are represented as coherent superposition of lp-lh states. In stable nuclei, 
because the Fermi energies are deep, lp-lh states between tightly bound single-particle states having 
similar spatial characters only contribute, and these lp-lh states concentrate around the nuclear surface. 
Around neutron drip line, by contrast, lp-lh states between tightly bound , loosely bound, resonance, 
and non-resonant continuum states, contribute. Because each single-particle state has different spatial 
extent, lp-lh states among them have different spatial characters. Therefore it is a non-trivial problem 
whether vibrational modes can be realized as a result of coherency between such lp-lh states. 

In Section|2]we examine spatial properties of lp-lh states in neutron drip line nuclei. First we study 
characters of lp-lh states without pairing correlations by solving the Woods-Saxon potential. Next, by 
solving simplified HFB equation in coordinate space, we emphasize that selfconsistent pairing correlations 
play important roles to realize collective motions in neutron drip line nuclei. In SecOHwe perform Skyrme- 
HFB plus QRPA calculations for the first 2 + states in neutron rich Ni isotopes. By comparing three types 
of calculations, HFB plus QRPA, resonant BCS plus QRPA, and RPA, we show that selfconsistent pairing 
correlations play important roles to realize low-frequency collective vibrational excitations in neutron drip 
line nuclei. 



2 lp-lh states in neutron drip line nuclei 
2.1 lp-lh states without pairing correlations 

We consider the spatial localization of lp-lh states as an important condition for collective vibrational 

motions. Before studying lp-lh states in neutron drip line region, we examine neutron lp-lh states in 

208 Pb as a typical example of stable nuclei for comparison. Because the Fermi energy is about -8 MeV, all 

wave functions around Fermi energy are tightly bound and spatially localized as shown in Fig^ These 

wave functions are provided by solving the Woods-Saxon potential. The remarkable point is that all 
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Figure 2: The spatial distributions of lp-lh states F^ h (r) for octupole transitions and the corresponding 
radii flg in 208 Pb are shown as a function of particle-hole energies s p — Eh- 

wave functions have large overlaps around the surface region. We define the spatial distributions of lp-lh 
states as 

Fih\r) = <p p (r)-r L - Vh (r), (1) 

where L is multipolarity of the transition and ip p (ifh) are radial wave functions of particle (hole) states. 
The spatial integration of (r) gives the radial part of the particle-hole transition matrix element 

(p|? ,L |/i) = J drFph\r). In FigEl F^\r) for octupole transitions are plotted. The hole (particle) 
states are indicated by lines (symbols), and each combination of a line and a symbol represents a lp- 
lh configuration. As seen, all lp-lh states are localized around the surface region. This feature is an 
important condition for realization of collective motions, and this condition is satisfied in stable nuclei. 
For further discussion, we define the "radii" of lp-lh states R^J as the second moment of F^'(r), 



2 



T 




Vws [MeV] 

Figure 3: The neutron single-particle energies around 86 Ni are plotted as a function of the depth Vws- 
Vws ~ -41 MeV corresponds to 86 Ni. 

In FigEJ R$ are plotted as a function of particle- hole energies e p — Eh- All concentrate around the 
surface region in stable nuclei. 

Next we investigate how the distributions of R p ^ change around neutron drip line region. We consider 
neutron rich Ni isotopes as a typical example. In the Hartree-Fock (HF) calculation with Skyrme SLy4 
parameter, 86 Ni is the neutron drip line nucleus. By taking into account pairing correlations, more 
neutrons can bound. However the precise position of the neutron drip line depends on the treatment of 
pairing correlations. For example, the drip line nucleus is 88 Ni in Ref.j3] and 92 Ni in Ref.j^. The predicted 
drip line also depends on the effective interactions and the frameworks, i.e., relativistic or non-relativistic 
approaches (for example, In the present study we consider up to 88 Ni for further discussion, because 

our purpose is to investigate the qualitative aspects of vibrational excitations in neutron drip line nuclei. 

In Fig|3 the neutron single-particle energies around 86 Ni are shown as a function of the depth Vws 
of the Woods-Saxon potential. Vws ~ — 41 MeV corresponds to the single-particle energies in 86 Ni 
calculated by HF with SLy4. We consider lp-lh states of quadrupole transitions around loosely bound 
3si/ 2 state (the Fermi level). The resonant states, 2d 3 / 2 and I57/2, are particle states. In Fig0]the radii 
of single-particle levels Ri are shown. The radii of bound states with large-£, l<7g/2 and 2d§/2, don't 
change as approaching the drip line. On the other hand, the radius of 3si/ 2 state becomes rapidly large 
in the limit of zero binding energy, i.e., neutron halo. Concerning particle states, 2c? 3 / 2 and I57/2, while 
the single-particle energies are negative, these radii change slowly as a function of Vws- On the other 
hand, once the single-particle energies become positive (resonant states), the radii increase suddenly. In 
Fig^the radii of lp-lh states R p h for quadrupole transitions are plotted as a function of Vws- The hole 
(particle) states are indicated by lines (symbols). Up to Vws < —44 MeV, corresponding to stable nuclei 
region, all R p h are localized around the surface region. By contrast, in the drip line region Vws > —44 
MeV, Rph strongly depend on the particle- hole configurations. Size of R p h is mainly determined by 
the spatial extent of hole wave functions. The localization of lp-lh states becomes more sparse and 
the condition for realization of collective motions is not satisfied as approaching the drip line. We may 
conclude that neutron drip line region is under the situation that collective motions are suppressed and 
single-particle like excitations appear dominantly. In next subsection we discuss how this situation is 
significantly modified by including selfconsistcnt pairing correlations. 

2.2 Roles of selfconsistent pairing correlations 

It is a well-known fact that pairing correlations are important ingredient to describe collective vibrational 
excitations. In this study we want to emphasize that selfconsistent pairing correlations described in the 
HFB theory play important roles in neutron drip line nuclei in non-trivial ways. 

Pairing anti-halo effect [S] is an important property of the HFB theory, and qualitatively modifies the 
asymptotic behavior of the quasiparticlc wave functions. In HF the asymptotic behavior of the bound 
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Figure 4: The radii of neutron single-particle levels Ri and the radii of lp-lh states R p h for quadrupole 
transitions around 86 Ni as a function of the depth Vws- 

states wave functions (ey, n < 0) for r — > oo is 

f?> F ( £ ij,n, r) -> exp (-aij >n r) , (3) 



where aij >n = \J 1 —2meij^ n /%. The decay constant aij. n can approach zero for ey^ — > 0. For positive 
energy states (e > 0), 

(e, r) Cy [cos (%) rj, (fir) - sin (<5y) rn* (fir)] (4) 



where /3 = y/2me/h. ji and n; are spherical Bessel and Neumann functions, and <Jy is the phase shift 
corresponding to the angular momentum (Ij). 

By taking into account pairing correlations, irrespective of selfconsistent or non-selfconsistent pairing, 
quasiparticle wave functions have two components, hole wave functions vij(r) and particle wave functions 
uijif). In HF plus BCS calculations both Uy(r) and uy(r) are just proportional to HF single-particle 
wave functions, 

v lj ( £lj , n ,r)=v l j, n ^(e lj , n ,r), (5) 

where uy )fJ (lty,n) are occupation (unoccupation) amplitudes. Therefore both hole and particle wave 
functions have the same asymptotic behavior for r — > co. Eqs.JHJ) are usually defined only for bound 
states. Recently the HF plus resonant BCS method was proposed [71IH]. In this approach resonant states 
are take into account for particle states. Applicability of the resonant BCS method for ground state 
properties is discussed in Ref.JH]. In HFB, by contrast, the asymptotic behavior of the quasiparticle wave 
functions is qualitatively different. We examine the HFB equation in coordinate space [31 1101 ITT] . 

hij{r)-X A(r) \f uij (E qp ,r) \ _ ^ / uij(E qp ,r)\ ^ 



A (r) -hij (r) + A J \ »y (E qp , r) J ^ \ v u (E t 

where hij (r) and Aij (r) are the mean field and the pairing field respectively. For the discrete region 
(E qp + A) < 0, both hole and particle wave functions decay exponentially at infinity, 



uij (Eij, n ,r) -> exp f-aj^r) 
vij (Eij, n ,r) -> exp (-/^rj 



(7) 



where cty n = \f— 2m (i5y n + X)/h and j3y n = y/2m (Eij n — \)/%. For the continuum region (E qp + X) > 
0, 



uij (E, r) -> Cy [cos (<Sy) rj; (a (2 V) - sin (dy) rn z (a (2) r)] 
t,y(£,r)^cxp(-/?( 2 )r) 4 



(8) 



where a (2) = ^2m (E + X)/h and (3^ = y/2m(E-X)/h. The hole wave functions always decay ex- 
ponentially. The minimum value of the decay constant (3 is estimated as following. If a discrete quasi- 
particle state with minimum energy is present, the minimum value of the decay constant becomes 

ffl = y/2m (E^-Xj/h. The HFB quasiparticle energies are well approximated by the canonical 

quasiparticle energies [H]. E r ~ E^ an = - A) 2 + (A™™) 2 , where e£°™ and A™ n are the canonical 

single-particle energy and the canonical pairing gap. This means that (3^ stays at finite value 



ffi _> yJImE^h, > ^2mA^ an /h > (9) 

with finite A° an for A — > 0. Selfconsistent pairing correlations affect the asymptotic behavior of all hole 
wave functions, and especially act against a development of an infinite rms radius that characterizes the 
wave functions of s and p states in the limit of vanishing binding energy. We call this localization of 
hole wave functions " pairing anti-halo effect" . Thought this terminology was used for the localization 
of the normal density distribution p(r) in Ref.|S], the microscopic mechanism is the same. In neutron 
drip line region where discrete HFB solutions don't exist, the analysis is more complicated. However, 
as numerically shown in Ref. 6 , pairing anti-halo effect plays the roles to localize hole wave functions in 
similar way. 

We solve the coordinate space HFB equation Eq. (JSJ with box boundary conditions Rb ox = 30 fm. For 
the HF potential the Woods-Saxon potential is used. The pairing potential is selfconsistently derived by 
the density dependent pairing interaction, 
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p(r) 



S(r-r'). (10) 



The parameters are V pa i r = —350 MeV fm -3 and p c = 0.16/2 fm -3 with the cut-off energy E cut = 50 
MeV. p c is taken to be the half of the saturation density, because Eq.® is solved only for neutrons. This 
pairing parameter set gives the averaged pairing gap A 1.5 MeV in the region —51 MeV < Vws < —39 
MeV. The average pairing gap is defined as the integral of the pairing field with the abnormal density 

X = J drp{r)\{f) 
I drp{r) 

In Fig|SJ the radii of hole wave functions around 86 Ni provided by HFB calculations are plotted. By 
comparing the radii of the single-particle states without selfconsistent pairing correlations (HF and HF 
plus resonant BCS in FigQJ, the spatial distributions of the hole wave functions are completely different. 
Due to pairing anti-halo effect, the radii in HFB spatially concentrate around the surface region. 

We consider how spatial characters of lp-lh states are modified by taking into account selfconsistent 
pairing correlations. First of all, by including pairing correlations, irrespective of selfconsistent or non- 
selfconsistent, the possible lp-lh configurations increase considerably, i.e., not only particle-hole but also 
particle-particle and hole-hole configurations can contribute. We define a spatial distribution of lp-lh 
state from a hole part vij(E) to a particle part ui>j>(E'), 



F [L) [ui'3'{E'),vi {E)](r) = u lir {E' ir )-r L -v l0 {E,r), (12) 



and the radius, 



/?Wr« (F>\ , (F ,,\ jdrr"FW[u ll3l {E'),v l0 {E)]{r) \ 1/2 

W '[ui>j>{E ),vij{E)\ = < — — — - > . (13) 

[ J dr FW[ui tjr (E'),vij(E)](r) J 

We use shorthand notations F^\r) and unless these don't cause confusions. In FigElthe radii 

R ( uv in HFB are plotted. As we have mentioned in subsec l2~Tl the radii of lp-lh states are mainly 
determined by the spatial extent of hole states. Because the radii R v in HFB are localized due to the 
pairing anti-halo effect, the radii R uv in HFB concentrate around the surface region. This localization 
of R uv can happen by treating pairing correlation selfconsistently. In resonant BCS, such localization 
can not appear. Generally speaking, the spatial distribution of particle-particle configurations between 
resonant states are always far outside from the surface region in resonant BCS. 
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Figure 5: The radii of hole wave functions R v and the radii of lp-lh states R v 
a function of the depth Vws- 



in HFB around 86 Ni as 



3 HFB plus QRPA calculations 

3.1 Formulation 

We consider the first 2 + states in neutron rich Ni isotopes by Skyrme-HFB plus selfconsistent QRPA 
calculations. By selfconsistent we mean that the HFB mean fields are determined selfconsistently from 
an effective force and the residual interaction of the QRPA problem is derived from the same force. The 
QRPA problem is solved by the response function method in coordinate space. A detailed account of 
the method can be found in Ref. |13ll*PIlll5| . Here, we just recall the main steps of the calculation. The 
QRPA Green's function G is a solution of a Bethe-Salpeter equation, 



G = G + G VG 



(14) 



The knowledge of G allows one to construct the response function of the system to a general external 
field, and the strength distribution of the transition operator corresponding to the chosen field is just 
proportional to the imaginary part of the response function. The residual interaction V between quasi- 
particlcs is derived from the Hamiltonian density (H) of Skyrme interaction by the so-called Landau 
procedure, 

Op /SO p & 

The index a (a = 1,2,3) stands for particle-hole (ph), particle-particle (pp), and hole-hole (hh) channels. 
The notation a means that whenever a is pp (hh) then a is hh (pp). The residual interaction V has an 
explicit momentum dependence, 



V(r, r') = F(Eu + *& v , ~&u + ~& v , ± *Vv, ± V v ]S(r - r') 



(16) 



The explicit form of the form factor F is shown in R.ef.|13j. These momentum dependence are explicitly 
treated in our calculation. Because we calculate only natural parity (non spin-flip) excitations, we drop 
the spin-spin part of the residual interaction. The Coulomb and spin-orbit residual interactions are also 
dropped. 

The ground states are given by Skyrme-HFB calculations. The HFB equation is diagonalized on a 
Skyrme-HF basis calculated in coordinate space with box boundary conditions Rbox = 20 fm. Spherical 
symmetry is imposed on quasiparticle wave functions. The quasiparticle cut-off energy is taken to be 
E C ut — 50 MeV, and the angular momentum cut-off is i m ax — 7h in our HFB and QRPA calculations. 
The Skyrme parameter SLy4 are used for the HF mean-field, and the density-dependent, zero-range 
pairing interaction Ea. (|10|l is adopted for the pairing field. The parameters are V pa i r = —555 MeV fm -3 , 
p c = 0.16 fm -3 , so as to give the average neutron pairing gap A„ s» 12/\/^4 in 86 Ni. 
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Figure 6: The B(E2, {f[ — > 2^) values and the excitation energies of the first 2 + states in neutron rich Ni 
isotopes calculated by HFB plus selfconsistent QRPA, resonant BCS plus QRPA, and RPA with Skyrme 
SLy4 force. 



3.2 Ni isotopes 

In order to emphasize the importance of pairing anti-halo effect for realization of collective vibrational 
excitations, we compare three types of calculations, HFB plus QRPA, resonant BCS plus QRPA (no 
pairing anti-halo effect), and RPA (no pairing) for the first 2 + states in neutron rich Ni isotopes. In FigHJ 
the B(E2,0^ — > 2+) values and excitation energies of the first 2 + states in neutron rich Ni isotopes up 
to 88 Ni in HFB plus QRAP, and 86 Ni for resonant BCS plus QRPA corresponding the neutron drip line 
are shown. To construct the quasiparticle wave functions without pairing anti-halo effect, we perform 
resonant BCS calculations. In resonant BCS calculations for the ground states, we include resonant states 
(localized states) for particle states. To identify the resonant states in calculations with box boundary 
conditions Rbox = 20 fm, we introduced an additional cut-off radius Rbcs = 10 fm and the single-particle 
orbits having the radii less than Rbcs are taken into account for the pairing problems of the ground 
states. The resonant l(?7/2 and 2<i 3 /2 states are included. The orbits with the radii larger than Rbcs are 
included only for QRPA model space. Because we intend qualitative discussion of microscopic aspects of 
vibrational excitations, we performed constant-A BCS calculations to avoid technical complexities that 
are inherent in non-selfconsistent calculations (lack of selfconsistency in the pairing potential, different 
energy cut-off scheme for HF plus resonant BCS calculations E cut = 5 MeV and QRPA calculations 
E C ut = 50 MeV, and neglecting the contribution of non-resonant continuum state that plays important 
roles in pairing correlations in neutron drip line nuclei (Hl ll7llT%| ). The constant pairing gap A = 0.8 MeV 
is chosen to reproduce low-lying quasiparticle energies in Skyrme-HFB calculations. We have checked the 
sensitivity of the pairing gap parameter, and the calculated excitation energies and the B{E2) values are 
qualitatively the same. For the residual interaction we used the exactly same interaction in HFB plus 
QRPA calculations. 

From 80 Ni to 82 Ni, HFB plus QRPA and resonant BCS plus QRPA calculations give qualitatively 
similar results. The energies decrease and the B(E2) values increase. On the other hand, as approaching 
the neutron drip line, the calculated B(E2) values exhibit different behavior. In HFB plus QRPA 
calculations the B{E2) values increase monotonically. On the other hand, the B(E2) values suddenly 
decrease from 82 Ni to 86 Ni in resonant BCS plus QRPA. The main qualitative difference between HFB 
plus QRPA and resonant BCS plus QRPA is the presence of pairing anti-halo effect. In resonant BCS, 
due to lack of pairing anti-halo effect, particle-particle configurations between resonant states, I57/2 and 
2<i 3 / 2 , and also particle-hole configurations with hole 3si/2 state are spatially decoupled to the other lp-lh 
states, and can not participate in the collective excitations. The results of HF plus RPA calculations for 
84 Ni and 86 Ni are also shown. These excitations have single-particle like characters, and the excitation 
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3/2 



energies are almost the same with specific lp-lh energies (2d 5 / 2 — > Ssi/2 m Ni and 3s\/2 
in 86 Ni). The B{E2) values in resonant BCS plus QRPA are larger than the values in HF plus RPA, 
because hole-hole configurations can contribute to make the collective modes in resonant BCS plus QRPA 
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calculations. 

For another example that indicates the importance of selfconsistcnt pairing correlations, we have 
shown that neutron pairing correlations play important roles to explain the experimentally observed 
anomalously large B(E2) values in neutron rich nuclei, 32 Mg and 30 Ne |13| . This is an actual example 
that pairing anti-halo effect plays important roles to realize collective motions in neutron drip line region. 

4 Conclusion 

We discussed the important microscopic aspects of low-frequency collective vibrational excitations in 
neutron drip line nuclei. We emphasized that pairing anti-halo effect in the HFB theory plays important 
roles to realize collective motions in loosely bound nuclei. We performed Skyrme-HFB plus selfconsistent 
QRPA calculations for the first 2 + states in neutron rich Ni isotopes. By comparing three types of 
calculations, HFB plus QRPA, resonant BCS plus QRPA, and RPA, the importance of pairing anti-halo 
effect was shown in realistic calculations. 
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